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We calculate the CP-averaged branching ratios and CP- violating asymmetries for the two- 
body charmless hadronic decays A;, — > pn, pK in the perturbative QCD (pQCD) approach 
to lowest order in as- The baryon distribution amplitudes involved in the factorization 
formulae are considered to the leading twist accuracy and the distribution amplitudes of the 
proton are expanded to the next-to- leading conformal spin (i.e., "P" -waves), the moments 
of which are determined from QCD sum rules. Our work shows that the contributions from 
the factorizable diagrams in A^ p-K, pK decays are much smaller compared to the non- 
factorizable diagrams in the conventional pQCD approach. We argue that this reflects the 
estimates of the p transition form factors in the kx factorization approach, which are 

found typically an order of magnitude smaller than those estimated in the light-cone sum 
rules and in the non-relativistic quark model. As an alternative, we adopt a hybrid pQCD 
approach, in which we compute the factorizable contributions with the —t p form factors 
taken from the light cone QCD sum rules. The non- factorizable diagrams are evaluated 
utilizing the conventional pQCD formalism which is free from the endpoint singularities. 
The predictions worked out here are confronted with the recently available data from the 
CDF collaboration on the branching ratios and the direct CP asymmetries for the decays 
Af, — !■ piT, and A^ —t pK . The asymmetry parameter a relevant for the anisotropic angular 
distribution of the emitted proton in the polarized A^ baryon decays is also calculated for 
the two decay modes. 

PACS numbers: 14.20.Mr, 12.38. Bx, 12.39. St, 13.30.Eg 



2 



I. INTRODUCTION 

The motivation to investigate b quark decays is attributed to their sensitivity to the quark 
flavor structure, which leads to an extremely rich phenomenology, studied mostly in the context of 
S-meson decays. However, heavy baryons containing a 6-quark have been observed at the Tevatron 
and they will be even more copiously produced at the large hadron collider (LHC). Their weak 
decays may provide important clues on the flavor-changing currents beyond the standard model 
(SM) in a complementary fashion to the B meson decays. A particular advantage of the bottom 
baryons over S-mesons is their spin, which provides a unique way to analyze the helicity structure of 
the effective Hamiltonian for the weak transition in the SM and beyond. Also, such baryon decays 
are flavor self-tagging processes which should make their experimental reconstructions easier. 

Theoretical analysis of non-leptonic decays are based on factorization theorems, which are the 
fundamental tools of the QCD perturbation theory enabling the separation of physics at different 
energy scales. The theoretical basis of the factorization theorem is a generalization of the Euclidean 
operator product expansion to the time-like domain. The proof of the factorization theorem has 
been worked out using the perturbative QCD approach based on the analysis of Feynman diagrams 
in the so called Collins-Soper-Sterman (CSS) formalism 

m 

. Equally importantly, the large 

mass of the heavy quark makes the formidable strong interactions effects controllable and they can 
be studied systematically using methods based on heavy quark expansion. 

The basic formula for the calculation of the branching ratios for the decays of the A;, baryon 
into two light hadrons is based on an operator realization of the diagrammatic analysis which can 
be described most easily for the calculation of the hadronic matrix element of B mesons decays 
into two light hadrons hi and /12. With the insertion of a set of the weak interaction operator Oj 
between the initial B meson and the final decay products hi and /12, the decay matrix element is 
obtained from the following formula {4] 

{hih2\Oi\B) = <^>h,{u) (T^(^x)F^'^nO) +C^^(r,n) H^'^i(r,0)) (1) 

involving the QCD form factor F^^^ (0) and an unknown, non-local form factor r.^^^ (r, 0) at the 
leading power in the A/m^ expansion. Different treatments of the various parts in the factorization 
formula ([T|) have led to three popular theoretical approaches to study the dynamics of non-leptonic 
two-body B meson decays, which are known as the perturbative QCD (pQCD) Q], QCD fac- 
torization (QCDF) and SCET approaches [3,1^) 9|. The function H^'*i(r, 0) is supposed to be 
dominated by perturbative hard-collinear interactions, and can be further factorized into light-cone 
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distribution amplitudes ^^l"^)) ^/ii(^) and a jet function J{t;lo,v) 

H^^i(r,0) = J{t;uj,v) $5(0;) <^hAv), (2) 

when the hard-colhnear scale ^JrnijKQCD is integrated out 

In contrast to these two latter approaches based on the collinear factorization theorem, the 
pQCD approach, which is developed in the framework of factorization, is free of the singularities 
from the end-point region of the parton momentum fractions. The pQCD approach has been widely 
applied for the calculation of the non-leptonic two-body B decays and it has proved itself to be 
successful in the description of exclusive processes with typical momentum-transfer of a few GeV. 
A hallmark of this approach is that the form factor F^^^ (0) is assumed to be dominated by short- 
distance contributions and it is therefore calculable in the perturbative theory. Soft contributions, 
though playing a role, are less important because of the suppression from the Sudakov mechanism 
embedded in the kx and threshold resummations . Current applications of the kx factorization 
theorem to exclusive processes are restricted to the leading order (LO) in the strong coupling 
constant a^. In this, the infrared divergences involved in the radiative corrections to the weak 
transition vertex are absorbed in the hadronic distribution amplitudes in a gauge invariant manner. 
The factorizable, non-factorizable and power-suppressed annihilation contributions are calculable 
in this framework free of the end-point singularities. 

In the case of non-leptonic two-body B decays, the decay matrix elements, in most cases, are 
dominated by the factorizable term, i.e. the first term on the right-hand side of Eq. ([T]), whereas 
the second term, the non-factorizable one, produces a perturbative correction. Since the first term 
proportional to the form factor F^^^ (0) is in pQCD very similar to the other approaches, where the 
form factors are input, the pQCD approach gives in most cases similar results for the non-leptonic 
B decays as the other two approaches mentioned above, though there are differences in detail. 

In the application of pQCD to two-body non-leptonic heavy baryon decays, we do not expect 
a similar pattern as in the non-leptonic B meson decays on general grounds. In particular, in the 
analysis of the hadronic decays of baryons, a large number of Feynman diagrams contribute to the 
hard amplitudes even in the lowest-order. Taking the — > pvr decay as an example, some 200 
Feynman diagrams need to be calculated as can be seen in section IIIIl These diagrams involve 
the exchange of two gluons, involving topologies where both gluons are attached to one of the 
light quarks emerging from the weak interaction vertex. As some of these diagrams build up the 
transition form factor, they receive contributions in a^, yielding small values for them. Another 
challenge for the baryonic transition is that the light-cone distribution amplitudes (LCDAs) of 
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the baryons are less known in the hterature. LCDAs are fundamental non-perturbative input 
to regularize the infrared divergence appearing in the radiative corrections in the factorization 
formalism of the pQCD approach. In view of this, applications of the pQCD approach to non- 
leptonic two-body 6-baryon decays are not worked out to a satisfactory level, and hence this area 
is essentially an uncharted territory. 

A first attempt to apply the pQCD approach to the baryonic transitions was made in 12l |. 



where the proton Dirac form factor is calculated taking into account the Sudakov suppression 
resulting from the resummation of the large double logarithms involved in the radiative corrections. 
Subsequently, the proton form factor was recalculated in [3] by refining the choice of the evolution 
scale of the proton wave functions and the infrared cutoffs for the Sudakov resummation, which lead 
to predictions for the Dirac form factors which are consistent with the experimental data. Following 
Refs. QQ^he se,„ilepto„ic charmless decays A. Q. the se.ileptonic cha™i„, decay 

Afe AJv [151 . Il6l |. the radiative decay A;, A7 [171], and the nonleptonic charming decay 
Aft — > KJ/if) (l8| have been investigated in the framework of the /cr factorization scheme. However, 
a study of the charmless hadronic decays Af, /ii/i2) which has been undertaken in the generalized 
factorization approach [3], to the best of our knowledge, is still lacking in pQCD. Our aim is to 
fill in this gap and provide further tests of the knp factorization formalism to gain insight on the QCD 
dynamics of these decays. In doing this, we have included the current information on the CKM 
matrix elements, updated some input hadronic parameters, such as the distribution amplitudes 



2l| making use of the conformal symmetry of 



of the proton, which are systematically studied in 
the QCD Lagrangian, and have used data to fix some other input quantities. We find that the 
non-factorizable contributions to the hard amplitudes overwhelm the ones from the factorizable 
diagrams in the baryonic decays A5 ^ pvr, pK. This feature of the 6-baryonic decays is at variance 
with what is found in the naive factorization approximation and in the corresponding two-body 
B meson decays. Large non-factorizable effects existing in the charmed baryon decays have been 
pointed out in the literature ii2^, where it is observed that the non-factorizable diagrams escaping 
from the helicity and color suppression can be comparable to and sometimes even dominate over 
the factorizable contributions. 

The layout of the paper is as follows: In section II, we briefly review the pQCD approach 
and give the essential input quantities that enter this approach, including the operator basis used 
subsequently and the LCDAs for the pseudoscalar mesons, the proton as well as the A;, baryon. 
Input values of the various mesonic decay constants and the baryonic wave function at the origin in 
configuration space are also collected there. Section III contains the calculation of the A;, pir, pK 
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decays, making explicit the contributions from the external W emission diagrams T, the internal 
W emission diagrams C, the W exchange diagram E, the bow-tie contraction diagrams B and 
the penguin diagrams P, as shown in Fig. [3l Details of the calculations are relegated to the 
two Appendices (Appendix A, where the Fourier integration to derive the hard amplitudes in the 
impact parameter (or b) space are displayed, and Appendix B, where the factorization formulae for 
the Feynman diagrams corresponding to various toplogies are given) . The decay amplitudes called 
fi and /2, defined in Eq. ()44p . resulting from the diagrams with different topologies evaluated 
in the conventional pQCD approach are given numerically in Table IIIII We find that the T 
diagrams dominate the A;, pn, pK decays, as expected. Numerical values of the factorizable 
and non-factorizable contributions from the T diagram amplitudes /j(Af, — > pn, pK);i = 1,2, in 
the conventional pQCD approach are given in Table HVl From the entries in this table we observe 
that the factorizable amplitudes in these decays are essentially two orders of magnitude smaller 
than the corresponding non-factorizable amplitudes. The form factor gi responsible for the A^ — > p 
transition evaluated in various theoretical approaches are collected in Table |Vl and we find that 
gi calculated in the pQCD approach is typically an order of magnitude smaller than in other 
approaches 



191, 



23( 1 ■ where the form factors are dominated by soft dynamics. Subsequently, we 
employ a hybrid prescription to deal with the hadronic A;, pn, pK decays. In this approach, 
the factorizable contributions are parametrized in the naive factorization approximation, and the 
variation of the renormalization scale is assumed to reflect the effect of the vertex corrections. The 
non-factorizable diagrams are evaluated, as in the conventional pQCD approach, in the framework 
of the kT factorization. Following this procedure and utilizing the from factors calculated in the 
light-cone sum rules (LCSR), we reanalyze these two channels and give the numerical results for 
the amplitudes /i(Ab pn, pK);i = 1,2, for the factorizable and non-factorizable contributions 
from the hybrid scheme in Table IVIII We note that the factorizable contributions are much larger 
in the hybrid scheme and they constitute a good fraction of the corresponding non-factoriazable 
amplitudes. Numerical results for the charge-conjugated averages of the decay branching ratios, 
direct CP-asymmetries and polarization asymmetry parameter a are_ tabulated in Table IVTnl A 
comparison of our predictions with the available experimental data 



2j] are also included in this 



table. Section IV contains our conclusion and an outlook. 
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II. CONVENTIONS, INPUTS AND SOME FORMULAE IN PQCD 



A. EfTective Hamiltonian 



We specify the weak effective Hamiltonian [251]: 
Gf ' 
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We// = -^|K6K;[c^i(/^)Qi(/«) + c2(/i)g^(/x)] - w,[ECi(/«)ft(/^) 

where q = d,s. The functions (« = 1, 10) are the local four-quark operators: 

• current-current (tree) operators 

Ql = {Uabi3)v-Aiqi3Ua)v-A, Q2 = iUaba)v-Aiqi3U(^)v-A, 

• QCD penguin operators 

Qs = {qaba)v-A^{q'0'p)v~A, Qi = {q(3ba)v-Aj2(^al'l3)v-A, 



+ h.c. 



Q5 = {qaba)v-Aj2(l'0P^'^+'^^ 
q' 

• electro-weak penguin operators 



Qe = iqpba)v-AY,(^a(l'l3)v+A, 

q' 



Ql = 7;{lo,ba)v-A^eq'{qfjq'f3)v+A, Qs = -{qf5ba)v-A^eq'{q^q'p 



Q9 = ^{qabc)v-A^eq,{q'pq'p)v^A, Q 



'V+A, 
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g' 



aW-Aj^^q'^l'c^'lP' 
q' 



V-A, 



(3) 



(4) 



(5) 
(6) 



(7) 
(8) 



where a and /? are the color indices and q' are the active quarks at the scale nif,, i.e. q' = (n, d, s, c, b). 
The left handed current is defined as {q'^q'p)v-A = q'alvi^ ~ l^)q'i3 ^-^d the right handed current 
^-s {q'aq'f3)v+A = q'a1u{^ + l^)q'f3- For later applications it will be convenient to use the following 
combinations of the Wilson coefficients Qi 26l |: 



ai = C2 + Ci/3, 03 = C3 + C4/3, as = C5 + Cg/S, aj = Cj + Cg/S, 09 = ^9 + Cio/3, 

a2 = Ci + C72/3, 04 = C4 + C3/3, ae = Ce + C5/3, ag = Cg + C7/3, aio = Cio + C9/3, (9) 

where the scale dependence for the Wilson coefficients has been suppressed here. For convenience, 
we have given the combinations of the Wilson coefficients at three different values of the energy 
scale in Table. HI 
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TABLE I: Numerical values of the effective Wilson coefficients defined in the text at three different scales 
fj,, where mi, is taken as 4.8 GeV. 



M (GeV) 


0.5m(, 


mb 


1.5m(, 




l.()(i 


1. ():-') 


1.02 


a2(xl()-^) 


0.40 


10.:-! 


14.8 


a:i(x !()-■') 


().41 




2.()3 


a4(xlO 


—32.6 


—22.8 


— 18.3 


a5(xl0-3) 


-5.87 


-2.29 


-1.20 


a6(xl0-3) 


-48.2 


-29.8 


-22.5 


a7(xl0-^) 


12.6 


12.2 


12.0 


a8(xlO-4) 


9.79 


7.57 


6.69 


a9(xl0-^) 


-84.5 


-82.2 


-81.4 


aio(xlO-4) 


-0.32 


-8.20 


-11.9 



B. Kinematics 

The kinematic variables of the initial and final hadrons can be defined as follows. The A{, baryon 
is assumed to be at rest, and the proton recoils in the minus z direction, p, p' and q = p — p' denote 
the momentum of the baryon, the proton and the light meson, respectively. The momenta of 
their valence quarks are parametrized as 

p=(p+,p-,0) = ^(1,1,0), 

fc^ = (xip+,jo",kiT) , A;2 = (x2P+,0,k2r) , ^3 = (a^sP^, 0, ksr) , 
p' = (0,p'-,O) = (0,p-,O) , 

k[ = (0, x'lp'-, k'lT) , k'2 = (0, x'2p'-,k'2T) , k'^ = (0, xy-, k'ar) , 
g = (g+,0,O) = (p+,0,O) , 

gi = (yg+,0,qT), (?2 = ((l-y)g+,0,-qT), (10) 

where ki {k'l) is the b (u) quark momentum, Xi (x'j) are their longitudinal momentum fractions, 
and k^^ are the corresponding transverse momenta, satisfying J2i k[^ = 0. y is the longitudinal 
momentum fraction carried by the quark in the emitted light meson and is its transverse 
momentum. The kinematics of the non-leptonic two body decays of A;, is described in Fig. 1. 
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k'. 



p' = ^(0,1,0) 



FIG. 1: Kinematics of the non-leptonic two-body decays of A;, in the pQCD approach. 



C. Distribution amplitudes of pseudoscaleir mesons 



The hght-cone distribution amphtudes for the pseudoscalar meson are given by [27|, |28| 



(P(P)|gXz)gi,(0)|0) 



V6 



dxe 



ixP-: 



(p^ix) 
6 



dxe 



ixP-: 



75 P<i)^{x) + 'y5mo(j)^{x) + mo75(?^ - 1)0^(2;) 



where 



^f,{x) 



^x{l-x)[l + 0A4Cl/'{t)], 



2^/6 



2^6 



[l + 0.43C2'/'(t)], 



^/6 

fK 



2^/6 



x(l - x)[l + 0.17C7f/^(t) + 0.115C2/^(t)], 
[l + 0.24C2'/'(t)], 



^[Cl''^{t) + {).?,bCl''^{t)], 



2^/6 



and the Gegenbauer polynomials are defined as: 



cl'\t) = t, 

cTit) = i(3t2 - 1), 



a/3 

a/3 ' 
(11) 



(12) 
(13) 
(14) 
(15) 
(16) 
(17) 



(18) 
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and t = 2x — l. The decay constants of these mesons are fixed as /tt = 130 MeV and fx = 160 MeV 
in our numerical calculations. 



D. Distribution amplitudes of baryons 



1. Distribution amplitudes of the Af, baryon 



The Lorentz structure of the baryon wave function can be simplified using the Bargmann- 
Wigner equation [29] in the heavy quark limit, where the spin and orbital degrees of freedom of 
the light quark system are decoupled. In the transverse momentum space, the wave function of 



the Afo baryon is defined as 



3d, 



3li 



1 [t1 <^^1 d^l ikrw,Jik 

2V2N, J l\ (2vr)3 
/At 



{0\T[bl,{0)uj,{w2)d';{w,)]\A,{p)) , 



(19) 



where b, u, and d are the quark fields, i, j, and k are the color indices, a, f3 and 7 are the spinor 
indices, C is the charge conjugation matrix, Ai,{p) is the A;, baryon spinor. The normalization 
constant corresponds to the value of the wave function at the origin in the configuration space. 
The numerical value /a^ = 4.28^Q'g4 x 10~^ GeV^ used by us is determined from the experimental 



data on the semileptonic decay A^ — > AJvi 



32( 1 . The quoted value (within the ztla range) is also 



in agreement with the ones estimated in the QCD sum rule (QCDSR) approach j33l . \34 

The phenomenological model for the distribution amplitude of the A;, baryon employed in this 
work is borrowed from 



35| 



■il;{xi,X2,xs) = NxiX2X3exTp 



Ml 



mt 



(20) 



2/32x1 2/32x2 2/32x3_ 

with the shape parameter /? = 1.0 it 0.2 GeV and the mass of the light degrees of freedom in the 
Afe baryon being mi = 0.3 GeV. The normalization 



[dx]llj{xi,X2,X3) = 1 , 



(21) 



ete set of three-quark distri- 



3g] in the heavy quark limit 



leads to the constant = 6.67 x 10^2 _ \Yg point out that the compL 
butions amplitudes of the A;, baryon has been investigated in Ref. 
and the renormalization-group equation governing the scale-dependence of the leading twist dis- 
tribution amplitude is also derived there. It is shown that the evolution equation for the leading 



twist distribution amplitude includes a piece associated with the Lange-Neubert kernel 



37t | which 
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generates a radiative tail extending to high energies, and a piece relevant to the Brodsky-Lepage 
kernel [33], which redistributes the momentum within the spectator diquark system. It is sufficient 



to limit the accuracy of the current pQCD analysis to the leading twist approximation due to the 
still large errors of the experimental data. 

The model for the twist-2 distribution amplitude for the A;, baryon proposed in 



361] is: 



-e"^ + a2C^^\2u - l)-e"t 



(22) 



with eo = 200t6o°MeV, ei = GSOlgooMeV and a2 = 0.333lo:333. In the above representation, uj is 
the total energy carried by the light quarks in the rest frame of baryon and the dimensionless 
parameter u describes the momentum fraction carried by the u quark in the diquark system. The 
normalization of tp^'^^ {uj , u) is 

"1 



ujduj I du4)^'^^{uj,u) = 1 



(23) 



For comparison, we translate Eq. (j20p in terms of the variables oj and u of Ref. 361 ]: 



Mi 



Nuj^u{l - u) 



ULO (1 — u)lJ 



Ma, 



Ma, 



xexp 



Ml 



mt 



mt 



Ma, Ma, ^ '^'H Ma, 



(24) 



The shapes of the LCDAs '0**^^°(w,u) and %l:^^^{uj,u), given in Eqs. (j22p and (j24p . respectively, 
are shown in Fig. [21 and the various curves show the dependence on the input parameters of the 
models. The variations of 02 in ijj^^^ {u:,u) does not play a significant role in the behavior of 
ip'^^^ [uj,u)^ since the second moment is suppressed by eo/ei, and so we have fixed 02 = 0.333. 
At this stage, it is difficult to select one or the other of these LCDAs. The harder spectrum of 
^CQM(^^^) in W (the sum of the energy of the two light quarks in the rest frame of the A^-baryon) 
also reflects in the inverse moments, which are more important for the dynamics. Following Ref. 



361 ] ■ we define the two inverse moments involving negative powers of the variables oj and u, the 



fractional quark momentum 



{{uju)-\Kuv)) ^ 



1 du 
u 



QCD/CQM 



{uj,u), {uj ^{Auv)) 



dco / dnV^Q^^/^«^(a 



where an additional energy cut u < Auv is introduced to guarantee that the moments are finite 
in the presence of a radiative tail. The values of ((wn)"^) and (i^^""^) for Auv = 2.5 GeV and 
// = 1 GeV are summarized in Table [HI We note from this table that the moments of the two 
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7 
6 

O 

^ 4 
's 

is 

o 

u 

8 2 
1 


0.5 1 1.5 2 

o) (GeV) 

FIG. 2: The functions i^^^^{uj,u) and m) plotted against w for the fixed value u = 0.5. The 

solid, dashed-dotted, dashed-double-dotted, dashed-triple-dotted, and dashed-quartic-dotted curves, peaking 
typically around oj = 0.8 GeV, describe the distribution amplitude ip'-^'^'^'^ {uj,u) with the values of the 
parameters (/3 = 1.0 GeV, to/ = 0.3 GeV), (/3 = 0.8 GeV, to/ = 0.30 GeV), (/3 = 1.2 GeV, to/ = 0.30 
GeV), (/? = 1.0 GeV, m, = 0.24 GeV), (/? = 1.0 GeV, to/ = 0.36 GeV), respectively. The curves peaking 
around oj — 0.4 GeV( red curves) correspond to the distribution amplitudes 4'^^^ {uj,u), where the solid, 
dashed-dotted, dashed-double-dotted, dashed-triple-dotted and dashed-quartic-dotted curves correspond to 
the values of the model parameters {eq = 0.20 GeV, ei = 0.65 GeV), (eq = 0.14 GeV, ei = 0.65 GeV), 
(eo = 0.33 GeV, ei = 0.65 GeV), (eo = 0.20 GeV, ei = 0.35 GeV), (eo = 0.20 GeV, ei = 1.30 GeV), 
respectively. 




TABLE II: Typical inverse moments defined in Eq. at a fixed energy cut off Ajjv = 2.5 GeV and 
= 1 GeV for the two LCD As i>'^^^{uj,u) and ij:^^'^'^ {ljj,u) discussed in the text. 





(Lj-i)[GeV-i] 


((tJu)-i)[GeV-i] 




^•""-0.67 
^•"'-0.15 


c qQ+2.22 
O.OO_2.07 

9 rq+0.42 



distribution amplitudes n) and ^*^^^(u;,n) are compatible with each other -within the 

errors on the model parameters (-which are large) , -with the central values of these moments shifted 
to lo-wer values for the il)'~^^^{u},u) LCDA. For the numerical calculations presented here -we use 
ip'~^^^{oj,u) -with the quoted errors on the model parameters. 
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2. Distribution amplitudes of the proton 



Similarly, the wave functions of the final state proton has in leading twist the following form 



N 



At 7 I 



8V2N, 



(26) 



ieeging next-to-leading conformal spin, one obtains the following twist-3 distribution amplitudes 



21 



4>^{Xi,fi) = 120x,X2X^[cl)l{n)+<t>+{fi){l-3x3) 

4)"^{Xi,^) = 120X1X2X3(X2 - Xi)(/)3 (;U) , 

(t)^iXi,n) = 120X1X2X3 (f)l{fl) - ^ (03 - 03 ) (m)(1 - 3X3) 



Here the moments of the distribution amplitudes for the proton are determined by 
with all the parameters fixed at the scale = 1 GeV as 



(27) 



1/ 



N\ 



(5.0 ± 0.5) X IQ-^GeV^ 



A'l = 0.38 ±0.15, 



Vf = 0.23 ±0.03. 



(28) 



It is easy to see that the above proton distribution amplitudes satisfy the following relations 



0^(xi,X2,X3) = 0^(X2,X1,X3) , 
(/)^(xi,X2,X3) = -0'^(X2,X1,X3) 
0'^(xi,X2,X3) = 0'^(X2,X1,X3) . 



(29) 



E. A brief review of the conventional pQCD approach 



Factorization of amplitudes is a fundamental tool of QCD perturbative theory to deal with 
processes involving different energy scales. Based on the kx factorization, the pQCD approach 
provides a framework which has been applied to hard exclusive processes. In this approach, hard 
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gluon(s) exchange is essential to ensure the applicabihty of the twist expansion, and soft contribu- 
tions are expected to be less important owing to the suppression by the Sudakov factor. This is 
the case for the transition form factors involving mesons. We would like to take the ki, ^ p transi- 
tion form factors as an example, first to illustrate the pQCD factorization theorem, and then offer 
quantitative estimates for this form factor to check if the soft contributions remain sub-dominant 
or not in the baryonic transitions. 

The factorization theorem states that the transition form factor can be expressed as the convo- 
lution of hadronic wave functions , tpp and the hard scattering amplitude Th 

F = J\dx][dx'] J [(fkr] J [d^h^r]Mx',kij.,p',ii)TH{x,x',MA^,kT,^T,l^)^J^,ix,kT,P,f^), (30) 

which is usually transformed to the impact parameter b space to perform the Sudakov resummation 
of the double logarithms involved in the radiative corrections to the hadronic wave functions 

F = J\dx][dx'] J [d%] J [d'h']Vp{x',h',p',fi)TH{x,x',MA,,h,h',fi)VA,ix,h,p,fi). (31) 

Here V\^{x,h,p, /j,) and Vp{x' ,h' ,p' , fi) are the Fourier transforms of the ipAf^{x,k.T,p, fi) and 
k^,p', /i), respectively. Radiative corrections to the hadronic wave function can generate 
a soft logarithm asln{Qb), whose overlap with the original collinear logarithm leads to a double 
logarithm agln^iQb). This type of large logarithm must be organized in order to ensure the validity 
of the perturbative expansion. Resummation techniques have been developed to deal with such 
double logarithms. The result is a Sudakov exponential exjp[—s{Q,b)], which decreases fast with 
increasing b and vanishes at 6 = 1/Aqcd- 

The expressions for the Sudakov evolution of the hadronic wave functions VAb{x,b,p, /j,) and 
Vp{x', h',p', jj) can be expressed as products of the Sudakov exponents s(6, Q) and reduced hadronic 
wave functions, denoted by 7'Aj(x,b,p, /i) and Vp{x' ,h' ,p' , fj,): 

3 



VAt,{x,h,p,n) = exp 
Vp{x',h',p',fi) = exp 



1=2 



where s{b, Q) is defined as 



^(1) 

s{b,Q) = -2^9^^ 



Ad) 



i=l 



-2a('-') + w(i-' 



VAt,{x,h,p,iJ,) , 
Pp{x',h',p',fi), 

^(2) 



(32) 



—1 



ln{2q) + 1 lii(25) + 1 



+ 



4/?f 4/?i 
\n^{2q) -ln^{2b) 



In 



In 



(33) 
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with 



g = ln[Q/(V2A)], 6 = ln[l/(6A)], 



(34) 



and the coefficients A^''-') and /3j are 



33-2nf 153 -19n/ 

Pi = T7i ' P2 = 



12 

^2 
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3 9 3 27 3 ^2 ^' 



(35) 



nj is the number of quark flavors and is the Euler constant. We will use the one- loop running 
coupling constant, i.e. we pick up the first four terms in the expression for the function s{Q, b). 

Apart from the double logarithms due to the inclusion of the transverse momentum, large single 
logarithms from ultraviolet divergences can also emerge in the radiative corrections to both the 
hadronic wave functions and the hard kernels, which are summed by the renormalization group 
(RG) method 



9 ^, .d 



Vp{x',h',p',ii) = -3-fqVp{x',h',p',i^), 



9 X 9 
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TH{x,x',MA„h,h',fi) = —-rqTH{x,x',MA„h,h',n). 



(36) 
(37) 
(38) 



Here the quark anomalous dimension in the axial gauge is 7^ = —as/n. In terms of the above equa- 
tions, we can get the RG evolution of the hadronic wave functions and hard scattering amplitude 
as 



VA^{x,h, p,fi) = exp 
Vp{x',h',p',iJ,) = exp 
TH{x,x',MA^,h,h',n) = exp 



KW fJ- 



8 ft' dfi ■ 

" J KW A* 



' Kw' fJ- 

17 /•* dfl 



X VAb{x,h, p,w), 
xVp{x',h',p',w'), 
xTH{x,x',MA„h,h',t) , (39) 



The factorization scales w and w' represent the inverse of a typical transverse distance among the 
three valence quarks of the initial and final states. The choices of w and w' are 



..111. , 111 



(40) 



with the variables bi and b'l defined as 
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with the other 6jS and defined by permutation. The introduction of the parameter n is done from 
the viewpoint of the resummation, since the scale kw^'\^ with k of order unity, is equivalent to w'^''^ 
within the accuracy of the next-to-leading logarithms 411]. The variation of k represents different 



partitions of the radiative corrections to the perturbative Sudakov factor and the non-perturbative 
wave function. The best fit to the experimental data of the proton form factor determines the 
parameter 1.14 Q. 

Furthermore, loop corrections for the weak vertex can also give rise to another type of double 
logarithms a^ln^Xj, which are usually factorized from the hard amplitude and resummed into the 
jet function St{xi) to smear the end-point singularity. It should be pointed out that the Sudakov 



42]. The following 



factor from threshold resummation is process- independent, and hence universal 
approximate parametrization is proposed in Ref. [43] for phenomenological applications 

with the parameter c ~ 0.3 determined from the best fit to the next-to-leading-logarithm thresh- 
old resummation in moment space. The threshold factor modifies the end-point behavior of the 
hadronic distribution amplitudes and forces them to vanish faster as x ^ 0. Collecting everything 
together, we arrive at the typical expression for the factorization formula of the form factor in the 
pQCD approach 

[dx] j[dx'] J[d^h] J[d^h']rp{x',h',p',w')TH{x,x',MA^,h,h',t)rA,{x,h,p,w)St{x^'^) 



X exp 



(43) 



Apart from the hard perturbative kernel Th{x, x' , ...), the same expression holds for the mesonic and 
baryonic transition form factors. As we shall see quantitatively below, the hard perturbative kernels 
entering the latter are parametrically suppressed compared to the former. Physical interpretation 
of the Sudakov factor is well known [44], namely it is a probability distribution function for emitting 
no soft gluons. When a quark is accelerated in QCD, infinitely many gluons are emitted. Hence, 
we may observe many hadrons (or jets) at the end if gluonic bremsstrahlung occurs. Therefore, 
the amplitude for an exclusive decay A;, to a light baryon and a light meson is proportional to the 
probability that no bremsstrahlung gluon is emitted. This is just the role of the Sudakov factor 
in the kx factorization. It is known that the Sudakov factor is large only for small transverse 
intervals between the quarks in the hadron. A large transverse interval implies that the quarks in 
the hadron are separated and hence less color shielded. Thus the Sudakov factor suppresses the 
long distance contributions for the decay amplitude. 
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III. CALCULATIONS OF BARYONIC DECAYS pvr, pK IN THE PQCD 

APPROACH 

Topological diagrams responsible for the decay of A;, to a light baryon and a light meson are 
presented in Fig. [3l In terms of the hard-scattering mechanism, the exchange of two hard gluons is 
needed to ensure that the light spectator quarks in the initial states turn out as collinear objects in 
the final state. With this, the various diagrams for the — > pir decays in the pQCD approach in 
the lowest order are displayed in Appendix B. Fig. IHshows the external W emission diagrams. Fig. 
[5] the internal W emission diagrams, Fig. [6] the W exchange diagrams, Fig. [7] the bow-tie diagrams 
and Fig. [8] the penguin diagrams. We also include diagrams containing the three-gluon-vertex 
displayed in Fig. [9l Their contribution is, however, about an order of magnitude smaller than that 
from the external W emission (T) diagrams, but it can be comparable to that of the internal W 
emission (C) diagrams. As for A;, — > pK decay, only Figs. HI EJ El and [9] contribute to the decay 
amplitude. 



A. General factorization formulae for A;, pn, pK decays 

The Af, — > pvr, pK decay amplitude M. is decomposed into two different structures with the 
corresponding coefficients fi and f2- 

M=p{p')[fi + f2j5]Mp)- (44) 

using the equation of motion for a free Dirac particle. Similar to the factorization formula for the 
form factors of the ^ p transition, the coefficients fi{i = 1,2) can be expressed as 

// = Gp-^fAj,^ J^ '^ f[Vx] f[Vby[as{tnfa\tn^^^^^ 

144 V3 m=v,A,T^ 

x^l\h,b',hq)eM-S^]. (45) 

Here, fl {i = 1,2) denotes the contribution to the coefficient fi by the "j't/i" diagram displayed 
in Fig. [31191 and are the corresponding Wilson coefficients. The hard function Vt^{h, b' ,bq) arises 
from the Fourier transformation of the denominators of the internal particle propagators in the jth 
diagram. The hard amplitudes Hj^^"'\x , x' , y) depend on the spin structures of the three valence 
quarks in the proton and the form factors fi^2- The integration measure involving the momentum 
fractions can be written as 

3 3 

[T>x] = [dx][dx']dy , [dx] = dxidx2dx35{l — ^Xj), [dx'] = dx'idx2dx'^5{l — ^x^), (46) 

1=1 i=l 
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u , u 

FIG. 3: Topological diagrams responsible for the decay Ab pn, where T denotes the external W emission 
diagram; C represents the internal W emission diagram; E labels W exchange diagram; B denotes the 
diagram that can be obtained from the E type diagram by exchanging the two identical down quarks in the 
final states; and P represents the diagram that can only be induced by the penguin operators. 



and the expressions for the measure of the transverse extent [Db] will be shown in the factorization 
formulae given in Appendix B. 

The exponents in the Sudakov factor are determined for the factorizable diagrams by 

S\x,x',b,b') = j2s{w,kt) + I f' '^^,{aM)+jls{w',k'r) + ^ , (47) 

j_2 J KW J KW' fJi 

and for the non-factorizable diagrams by 

S^{x,x',y,b,b',bq) = ^ A;+) + | /" ^-fq{as{fl)) + '^s{w' ,k'~) + 3 f ^jq{as{fl)) 

+ ^.K,g+) + 2 f'^^^iasm , (48) 
where P is the typical energy scale of the "jth" diagram and is chosen as 

=max{ti,t2,t-'^,t-l,w,w',Wg), (49) 
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TABLE III: The coefficients /i and /2 contributed by the Feynman diagrams with definite topologies in the 
A{, —^ pi: decay based on the conventional pQCD approach. 



/l /2 



T 


-2.42 X 10-9 - 


- i2.07 X 10 


-9 


-1.74 X 10-3 - 


- il.22 X 10 


-9 


C 


2.05 X IQ-io - 


i4.60 X 10- 


10 


-2.35 X 10-1° - 


hi4.77x 10 


-10 


E 


2.89 X 10-" - 


i8.95 X 10- 


12 


1.11 X 10-" - 


i4.36 X 10- 


12 


B 


-7.00 X 10-" - 


Ki3.33 X 10 


-10 


2.21 X 10-1° - 


i4.04 X 10- 


11 


P 


-6.84 X 10-12 - 


Fi4.85 X 10 


-11 


7.00 X 10-12 - 


i4.75 X 10- 


11 


G 


1.37 X 10-1" + 


il.71 X 10" 


11 


-1.60 X 10-1° - 


hi2.01 X 10 


-10 



where the hard scales t{, t2 are relevant to the two virtual quarks, and ti^, t\ are associated with 
the two hard gluons. w and w' have been given in Eq. (j40p and Wq = 1/bq. The maximum in the 
above choice simply indicates that the hard scales should be larger than the factorization scales. 

The factorization formulae for some typical diagrams corresponding to different topologies in the 
Ab — > pir decay are given in Appendix B. The corresponding factorization formulae for A^ pK 
decay can be obtained directly following the same rules. 

B. Numerical results for Af, ^ pir, pK decays 

For the CKM matrix elements, we use as input the updated results from j3] and drop the 
(small) errors on Vud-, Vus and Vtb'- 

\Vud\ = 0.974, \Vus\ = 0.225, |Kfe| = (S.SOl";!^) x IQ-^, 

\Vtd\ = (8.59t°;i) X 10-3, \Vts\ = (40.411?;?^) x IQ-^, \Vtb\ = 0.999, (50) 

/?=(21.58t°|l)°, 7 = (67.8tt^)°. 

It will be shown that the CKM factors mostly yield an overall factor for the branching ratios 
and do not introduce large uncertainties to the numerical results. 

We start by discussing the numerical results in the conventional pQCD approach. To that 
end, we list the coefficients /i and /2 defined in Eq. (j44p contributed by the Feynman diagrams 
with different topologies in the A;, — > pir decay in Table lllli From this table, we observe that the 
amplitudes satisfy the relations T ^ C ^ E. 

As mentioned earlier, the T type diagrams dominate the A^ — > pn decays. For this case we 
present the factorizable and non-factorizable contributions in the Af, — > pir decays in Table IIV[ 
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TABLE IV: The coefficients /i and /2 in the A;, pn, pK decays from the factorizabie and non-factorizable 
external W emission (T) diagrams in the conventional pQCD approach. 

factorizabie non-factorizable 
/i(Ah->p7r) 1.47X 10"" -il.97x 10"" -2.43 x 10"^ - i2.05 x 10"^ 
/2(Ah^p7r) 1.26X 10""-il.94x 10"" -1.75 x 10"^ - il.20 x 10"^ 
/i(Ab ^ pK) -1.52 X 10"" - iO.62 x 10"" -0.88 x 10"^ + iO.54 x 10"i° 
f2{Kb^pK) 0.17X 10"" -iO.eOx 10"" -1.06 x 10"^ + il.67 x 10"^ 



We observe that the factorizabie contribution is approximately two orders of magnitude smaller 
than the non-factorizable contribution. This is also the reason that the the conventional pQCD 



predictions for the semileptonic decay Af, plT> [39] and the radiative decay A^ — > A7 |17|] are 
much smaller than those evaluated in other theoretical frameworks (such as the constituent quark 
model or the QCD sum rules). 

The suppression of the factorizabie contributions in the conventional pQCD approach has been 
observed also in the analysis of the A;, KJ/ij: decays [31], where the non-factorizable contributions 
are also found almost two orders of magnitude larger than those from the factorizabie diagrams. 
In order to understand the large contribution of the non-factorizable diagrams in A;, decays, it is 
necessary to recall the role of the Sudakov factor in the kx factorization approach. As stated in 
section II, the Sudakov factor can only suppress the region with large Vs corresponding to small 
/ct's, and has almost no effect in the region where the transverse momentum Ut is large. Taking 
the non-factorizable diagram T25 as an example, the two virtual quarks can be on the mass shell 
even in the region with large kx- Therefore, this diagram is not subjected to the suppression from 
the Sudakov factor. It is then expected that the amplitudes for the non-factorizable diagrams 
should be much larger than those from the factorizabie diagrams, where the two virtual quarks 
can be on the mass shell only in the small kx region. Actually, a similar case also occurs in the 
hadronic B meson decays. There, the annihilation diagrams contributing to the B M1M2 decays 
in the pQCD approach is very important, which is responsible for the large CP violation and the 
enhancement of the transverse polarization fractions predicted in the kx factorization. The large 
contribution from the annihilation diagrams in the pQCD approach is due to the fact that the 
inner quark can be on the mass shell in the region of large kx- The numerical analysis also shows 
that the six non-factorizable diagrams Tig, T20, 721, T25, T31, T32 play the most significant role in 
the decay amplitude for the A;, pn transition. 
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TABLE V: The form factor 171 responsible for the At ^ p transition at zero momentum transfer, calculated 
by us (this work) and in the non-relativistic quark model (NRQM), LCSR, and in another pQCD approach. 
The uncertainties from the variations of the hard scale, Aqcd and the shape parameter /3 in the wave 
functions have been combined together in our work. 

NRQM \19\ LCSR (fuU QCD)[23] pQCD \l£ pQCD (this work) 
gi 0.043 0.018 2.3 x lO'^ 2.2t° j x lO'^ 



We consider the smallness of the factorizable contributions in the conventional pQCD approach 
as unrealistic. Consequently, we argue that the Ah ^ p transition form factors can not be reliably 
calculated in the perturbative scheme, i.e. these form factors are dominated by non-perturbative 
soft contributions, which can not be estimated in the pQCD approach. Of course, this could 
easily be checked by measuring the semileptonic A;, decays Af, pii^e, which depend only on the 
factrorizable diagrams. Pending this determination, we consider it as a more reasonable approach 
to calculate the Af, ^ p transition form factor by means of some non-perturbative method. 

The form factors Af, ^ p transition are defined as 

{p{p')\u7^J,b\Ah{p)) = p{p'){gii^,+ 92i(ypivq'' + g?,qpi)Ah{p), (51) 

where all the form factors gi are functions of the square of momentum transfer q^. We show in 
Table IVl numerical values for the vector transition form factor gi for the A^ — > p transition. These 
results are obtained in the non-relativistic quark model (NRQM) [l^, LCSR j^S], an earlier pQCD 
calculation [l4], and this work (also a pQCD calculation) for comparison. From this Table we 
see that the predictions for the transition form factor gi are scattered, with the NRQM [l^ and 



the LCSR [2^ values differing by a factor 2, but the two conventional pQCD results shown, while 
consistent with each other, are smaller from those obtained using the non-perturbative methods 
typically by an order of magnitude. 

To understand the marked difference of the form factor gi predicted in the pQCD approach 
and in the other frameworks, we recall that the hard dynamics is assumed to be dominant in 
the heavy-to-light transition form factors in the former and the soft contribution, which is not 
calculable, is assumed to be less important due to the Sudakov resummation. Table |V] suggests 
that the soft dynamics in the heavy-to- light transition form factors is the dominant effect, in all 
likelihood overwhelming the mechanism of the hard gluon exchange for the baryonic transitions. 
Similar large soft contributions have been also observed in the nonleptonic charmed meson decays 



21 



4g] as well as in the semileptonic — > A7, Al~^l~ decays [171 . |33|. It is found in [33|] that the 
hard contributions to the A;, ^ A form factors are almost an order of magnitude smaller than that 
those from the soft contributions. 

In the modified version of the pQCD approach, which we call hybrid pQCD, the form factors are 
taken as external inputs. The perturbative correction to the factorizable amplitude will then enter 
through the Wilson coefficients, which are known in next-to-next-to- leading order (NNLO), and 

)een recently calculated for the tree diagrams in the charmless 



the vertex corrections, which have_ 
hadronic B decays in NNLO 
guin amplitudes, are still not yet at hand, we follow the approximate (and less precise) approach 



48l |. As the complete NNLO corrections, including the QCD pen- 



proposed in Ref. [49(] to neglect the vertex corrections and vary the renormalization scale of the 
Wilson coefficients between O.Sm;, and LSm^. Surely, this step of the calculation can be system- 
atically improved once the complete NNLO virtual corrections are available. The non-factorizable 
contributions will be evaluated as already discussed in the conventional pQCD approach. 

Following the above procedure, we write the complete decay amplitude for A^ — > pir, pK as 

M{Ab^pTT) = Mf{Ab ^ ptt) + Mnfi^b ^ ptt) 

M {Ab ^pK) = Mf (Ab ^ pK) +Mnf [Ab ^ pK) , (52) 

where A4nf{Ab pir) and Mnf{Ab pK) denote the contributions from the non-factorizable 
diagrams and have been computed in the conventional pQCD approach. To calculate the factoriz- 
able amplitudes Aif{Ab pir) and Aif{Ab pK), we first need to deal with the hadronic matrix 
elements with the insertion of {V — A) {V + A) operators, i.e., the O5 — Os penguin operators. 
Making use of the Fierz identity, the factorization assumption and the Dirac equation, the matrix 
element of the operator Og can be written as 

:>M /i-/^, 11 A \ I vM 



{pM\Oe\Ab) = [Rf {p\q%b\Ab) + R^' {p\q%j5b\Ab)]{M\qj^{l - ^,)q'\0) (53) 



with 



rM = ?!!!M ^, rM ^ ?!!^M ^, (54) 

[mb - mu){mu + ruq) [mb + mu)[mu + mg) 

where the quark masses are the current quark masses. In addition to the form factors defined in 
Eq. (jSip . we need the matrix element describing the Ab ^ p transition induced by the axial- vector 
current 



(A(P) 1 57^756! Afe(P + q)) = A(P)(Gi7^ + G2ia^,q'' + Gsq^^MP + q). (55) 
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TABLE VI: Numerical values of the form factors gi and mt^^g2 at zero momentum transfer, responsible for 
the Ah ^ p transition, estimated in the LCSR approach [2^. 



form factors 


91 






0.018 


-0.159 



It is then straightforward to write the factorizable amplitudes A4/(Af, pn) and Aif{Ai) — > pK) 



as 



Mf{Ab^pTT) 

Gf 



V2 



fMp'] 



VubV*dai - VtbVid{a4^ + aio + i?i (ae + as)) 



5i(m2)(MA, - Mp) + g:i{ml)ml 



+ 



KfeKdOi - VtbVtdi'^i + aio - -^2(^6 + as)) 



Gi(m2)(MA,+Mp)-G3(m2)m2 



75 \Kb{p), (56) 



Mf{Kb^pK) 



VubV*gai - VtbV*g{a4 + aio + i?f (ae + as)) 



(7i(m|^)(MA, - Mp)+gsiml)ml 



+ 



- ^tbVtliai + aio - (ae + as)) 



Gi(m2,)(MA, + Mp) - G3{ml)m 



K 



75 Afe(p). (57) 



The masses of the pseudoscalar mesons of vr and K can safely be neglected, therefore only the form 
factors at the zero- momentum transfer will be involved in the numerical computations. 

To evaluate the A^ pir, pK decays numerically, we need to specify the form factors responsible 



for the Kb ^ p transition. As can be seen from Eqs. (j56ll57p . the form factors g^ and G3, whose 
contributions are proportional to the mass of the corresponding meson, are inessential for the 
calculation of the decay amplitudes. In view of the minor effects of these two form factors, it is 
quite adequate to determine them in terms of the relations derived in the heavy quark limit. As is 
well known, the form factors gi and Gi satisfy 



gi 
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Lri — 4i H 42, 



Go 



33 = (^3 = 



A. 



(58) 
(59) 



in the heavy quark effective theory (HQET), where the two independent form factors ^1 and .^2 are 
defined as 



(A(P)|6r5|A,(P + Q)) = A(P)[ei(g2)+^e2('z')]rA,(P + 



(60) 



with r being an arbitrary Lorentz structure and being the four-velocity of the Af, baryon. An 



analysis of the form factors gi and Gi has been performed in the LCSR 



23], which we shall use 
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TABLE VII: The coefficients /i and /2 in the decays A;, ^ p7r, pK contributed by the factorizable and 
non-factorizable external W emission (T) diagrams in the hybrid pQCD scheme. 

factorizable non-factorizable 

/i(Ah->p7r) 2.43 X 10^1° -i4.39 X 10^1° -2.43 x 10"^ - i2.05 x 10"^ 

/2(Ah^p7r) 2.64 X 10^1"- i6. 54 X 10^1" -1.75 x 10"^ - il.20 x 10^^ 

/i(Ab pK) -3.17 X 10-1° - il.22 x lO^i" -0.88 x 10"^ + iO.54 x 

f2{Kb^vK) 1.74 X 10-1" -il.96 X 10-1" -1.06x10-^ + ^1.67x10-9 



here. The numerical results for gi and m\^g2 needed for our numerical calculations are grouped in 
Table IVTl which correspond to = 0.050 and ^2 = —0.16. 

Utilizing the Wilson coefficients, the input form factors just discussed and the CKM factors 
given earlier, we can now compute the factorizable contributions to /i and /2 in the hybrid pQCD 
approach and compare them to the corresponding non-factorizable contributions, which have been 
given already earlier. The results are given in Table IVTIl From this table we see that the factorizable 
contributions are now much larger than in the conventional pQCD approach, though they are still 
smaller than the corresponding non-factorizable contributions. 

We are now in a position to present our final results concerning the branching ratios, direct 
CP asymmetries and the polarization asymmetry parameter a for the two decay channels in the 
conventional pQCD and in the hybrid pQCD approach. The CP-asymmetry ^cp(A|^ — > pT^~) is 
defined as follows: 

with Acp(AO pK^) defined similarly. The asymmetry parameter a associated with the 
anisotropic angular distribution of the proton emitted in the polarized Af, baryon decays is de- 
fined as follows: 

r = ro(l + ap-SAj (62) 
with p, sA|3 being the three-momentum and spin vector of the proton in the rest frame of the A^ 



baryon. The explicit expression of a can be written as 



^ 2kRe{ftf2) , . 

(l/iP + zi^l/sP)' ^ ^ 



with K = |p|/(-E'p + mp) = J {Ep + mp)/{Ep - nip). 
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TABLE VIII: The CP-averaged branching ratios, direct CP asymmetries and the polarization asymmetry 
parameter a for the Af, —^ pn, pK decays obtained in the conventional and the hybrid pQCD approaches. The 
errors for these entries correspond to the uncertainties in the input hadronic quantities, the scale-dependence, 
and the CKM matrix elements, respectively. Current experimental measurements at the Tevatron 
also listed. 



24j are 



pQCD (conventional) pQCD (hybrid scheme) 



Exp. 



BiAt- 

^cp(Ab 
Acp{Ab - 

a{Ab - 

a{Ab 



pit) 
pK) 
-^pn) 
^pK) 
PTt) 
pK) 



A f;fi+2-08-l-0. 70-1-0. 35 -.r.-^ 

74-0. 35-0. 35 ^ 

1 00-1-0.74-1-0.62-1-0.07 ^ i n- 

^•°^-0. 71-0. 80-0. 05 ^ 



-0.32 



+0.27+0.41+0.01 
-0.00-0.00-0.01 



r 91+2.42+0.30+0.42 -, r,- 
^■^^-1.89-0.10-0.37 ^ 

9 n9+0.78+0.55+0.10 i r,- 
^■'-'^_0.86— ^ Qn_n nc; ^ 

-0.31^ 



r, rvq + 0. 21 + 0. 13 + 0. 00 

"^■^"-•-0.00-0.04-0.00 



-0.83 



+0.03+0.03+0.01 
-0.01-0.07-0.01 



n nQ-i-0. 00+0. 00+0. 03 

^■""'-0.36-0.07-0.05 



-0.86-0.90-0.05 

1+0.28+0.32+0.01 



-0.00-0.00-0.01 

_n nc;+o-26+o.o3+o.oi 

U.UO_g^00-0.05-0.00 



-0.84 



+0.03+0.00+0.01 
-0.00-0.00-0.01 



"■'-'"-0.38-0.42-0.04 



3.5 ±0.6 ±0.9 X 10^6 

5.6 ± 0.8 ± 1.5 X 10^6 
-0.03 ±0.17 ±0.05 
-0.37 ±0.17 ±0.03 



We present our results in the two pQCD approaches and compare them with the current ex- 



perimental data from the Tevatron 



2j] in Table Ivnll The first error in the pQCD-based entries 



arises from the input hadronic parameters, which is dominated by the errors on the normahzation 
constant of the A;, baryon (taken as /a^ = 4.28lgg4 ^ 10~^GeV^) and the A^ baryon wave function 
shape parameter (taken as, (3 = 1.0 ± 0.2 GeV). The second error is the combined error from the 
hard scale t, defined in Eq. (j49p . which is varied from 0.75t to 1.25t, and the renormalization scale 
of the Wilson coefficients, given in Table I. The third error is the combined uncertainty due to the 
CKM matrix elements. 

We observe from Table IVlnl that the results for the conventional pQCD and the hybrid pQCD 
approaches do not differ very much, although in the hybrid approach the factorizable contributions 
have increased by almost an order of magnitude as compared to the conventional pQCD approach. 
The reason for this is that in the hybrid approach the factorizable amplitudes fi are still only a 
fraction of the non- factorizable amplitudes, as is apparent by comparing the results in Table IVTll 
Of course, it remains to be checked if the non-factorizable amplitude is correctly estimated in the 
pQCD approach for the 6-baryonic decays due to the exchange of two gluons. This involves, among 
other diagrams, those where both the gluons are attached to the same outgoing quark line (see, for 
example, the diagrams in the fourth row in Fig. 2]). These contributions are more sharply peaked, 
compared to the others encountered here or in the decays of S-mesons, which involve single gluon 
attachments on a quark line. 

The ratio of the decay rates for the A^ — > pvr and A;, pK decays, called below R^^KiAb), can 
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be calculated from Table IVIIIl and is estimated by us as 

in the hybrid pQCD approach. This can be understood from Eqs. (I56ti57p . which show that the 
QCD penguin operators contribute to the coefficients /i and /2 (defined in Eq. (I44p ) in the 
combination + RiQq and 04 — -R^^ag, respectively. This is quite different from the two-body 
hadronic decays of the B mesons, B — > PP or B ^ PV, where P{V) is a light pseudoscalar 
(vector) meson. The key point is that both the hadronic matrix elements (A(P)|s7^6|A;,(P-|-g)) and 
(A(P)|s7^756|Aft(P-|-g)) contribute to the baryonic decays. Theoretical predictions presented here 
deviate from the experimental data i?7r_R'(Ab) = 0.66 it 0.14 it 0.08 Whether this discrepancy 
reflects the inadequacy of the current theoretical formalism embedded in the standard model, or 
the standard model itself, or requires improved data remains to be seen. We note en passant that 
the estimates of the branching ratios for the decays A^ pvr and A^ pK, and hence of the 
quantity Rj^Ki-^b), reported in [19] in the generalized factorization approximation, are in error 
due to the incorrect relative sign of the two terms in Eq. (18) in that paper. We are convinced 
that the correct relative sign in question is given in our Eq. (j53p . 

As for the direct CP asymmetries, theoretical predictions suffer from large uncertainties due to 
the hadronic distributions, the hard scattering and the renormalization scales in the factorizable 
amplitudes. For the CP asymmetries, one needs the complete NNLO vertex corrections, as only 
with this input will it be possible to make quantitative predictions. As can be seen from Table 
IVIIH theoretical estimates for the parameter a for the decay Af, pn have negative values in both 
the pQCD approaches, reflecting the {V — A) structure of the weak current [5lJ- It is pointed out in 
52| that the parameter a in Bi{^~^) Bf{^^)P{V) decays approaches —1 in the soft pseudoscalar 



meson or vector meson limit, i.e., for mp ^ or my 0. This argument, however, is only valid 
for the tree-dominated processes. As for the A;, — > pK decay, the contributions from the QCD 
penguin operators are comparable to that of the tree amplitude. The operator Oq contributes to 
the Af, — > p transition via the {V + A) current (see Eq. (j53p ) and the Wilson coefficient og is very 
sensitive to the energy scale as can be seen from Table HI Hence, the asymmetry parameter a can 
flip its sign for the A^ — > pK decay due to the large penguin contributions. As a flnal remark, 
we find that the predictions for the parameter a in the A^ — > pn decay are relatively stable with 
respect to the variations of hadronic parameters, the CKM matrix elements and the hard scale, 



and therefore it serves as a good quantity to test the standard model 



3ii. 
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IV. DISCUSSIONS AND CONCLUSIONS 

Thanks to the current and impending experimental programs at the Tevatron and the LHC, 
dedicated studies of the decays of the A;, baryon (and other heavy baryons) will be carried out, 
following the first measurements of the decays A^ — > pvr, pK, performed at the Tevatron. Baryonic 
decays are flavor self-tagging processes. Therefore, they should be easier to reconstruct experimen- 
tally. In particular, the CP-asymmetry measurements amount to counting these self-tagged decay 
modes and their CP-conjugates. From the theoretical viewpoint, however, 6-baryon decays are less 
tractable as the underlying QCD dynamics is more involved. Hence, it is far from being obvious 
if the theoretical approaches developed for the quantitative studies of the two-body non-leptonic 
decays of the i?-mesons will work also for the corresponding 5-baryon decays. We have carried out 
an exploratory study of the charmless hadronic decays A^ pir, pK in the pQCD approach and 
find that the factorizable diagrams in the conventional pQCD approach contribute very little to the 
branching ratios, as the hard (perturbative) contributions to the baryonic transition form factors 
in this case turn out to be quite small compared to the estimates dominated by the soft dynamics. 
As an alternative, we adopted a hybrid approach, in spirit similar to the one advocated in Ref. [4^ 
for the analysis of the color-suppressed decays, such as J /iJjK^ . An essential characteristic 

of this hybrid scheme is that the transition form factors are treated as non-perturbative objects, 
i.e., they are input in the theoretical analysis and are not computed perturbatively, as in the con- 
ventional pQCD approach. Employing the form factors estimated in the LCSR approach, we find 
that the factorizable contributions are no longer negligible, though for the two decays worked out 
here, the amplitudes are still dominated by the non-factorizing contributions. 

Our predictions for the branching fraction for the decay Af, pir, which is dominated by 
the tree diagrams, are essentially in agreement with the current data, whereas estimates of the 
branching ratio for the A;, — > pK decay, dominated by the penguin-amplitude, are found to be 
smaller typically by a factor 2. This deserves an improved theoretical analysis, as the data gets 
consolidated. The asymmetry parameter a associated with the anisotropic angular distribution of 
the proton produced in the polarized A5 baryon decays is also derived and is found to be relatively 
stable with respect to variations of hadronic inputs and higher-order corrections in Af, pn decay. 
The asymmetry parameter a in the Af, pK decay, however, can flip its sign due to the large 
penguin contributions and the sensitive scale dependence of the effective Wilson coefficient ag (/x) . 
The Feynman diagrams (G) with the three-gluon-vertices present in the perturbative amplitudes 
included in this work are found to be less important compared with the T diagrams. However, 
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these three-gluon-vertex diagrams are comparable to the C diagrams, as can be seen from Table 
IIIH and hence they may induce significant corrections to the color suppressed modes, such as the 
Afc — > KJ/ij: decay. Finally, quantitative estimates of the CP asymmetries presented here show 
large scale uncertainties and require NNLO vertex corrections to be firmed up, which are not yet 
available completely. 
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APPENDIX A: FOURIER INTEGRATIONS AND fe-SPACE MEASURES 

We list below the Fourier integration formulae which have been employed in the derivation of 
the hard amplitudes in the impact parameter (or b) space. The symbols Ji, A'^i, Kq and Ki are 
the various Bessel functions; Zi are the Feynman parameters; and the relation 

Kn{-iz) = y e"^ [Jn{z) + iNn{z)] (Al) 

has been used used in the derivation of the Fourier transformation. With this, we get: 



(A2) 



(fc2 + I)'jl2 + jj) = ^ dz-^^{K,{^ib)9{Z,) + ^[N,{^\b) - .Ji(y^6)]0(-Zi)}(A3) 



/ 



oi(ki-bi+k2-b2) 



(kl + A){k^ + B)[iki + k2)^ + C] 



TT 



dzidz2 



TT 



Ki{VX^)9{Z2) + - 



lo Zi{l-Zl 

where vl > , and B, C are arbitrary , 



NiU X2\Z2\) - iJl{\ X2\Z2 



-Z2) , 



(A4) 
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/ 



gi(ki-bi+k2-b2+k3-b3) 



TT 



{kl + A){kl + B){kl + C)[{ki + k2 + k^Y + D] 
dzidz2dz3 



TT 



lo Zi{l - Zi)z2{l - Z2) \J\Z-i\ 

A,B > , and C, D arbitrary , 
with the variables, 



Ki(V^^)0(Z3) + - 



(A5) 



Z2 
X2 
Z3 
X3 



Az + B {1- z) , 

A{l-Z2) + 



Z2 



Zl(l - Zi 



■[B {l-zi) + Cz,] , 



/, , \2 -21(1 — -Zl) ,2 

{bi - zih2Y + — -h^ 



Z2 



2 ' 



Z2ZZ 



(A6) 



(A7) 



(A8) 



APPENDIX B: FACTORIZATION FORMULAE FOR THE FEYNMAN DIAGRAMS 

WITH VARIOUS TOPOLOGIES 



In this appendix, we would like to collect the factorizable formulae for typical diagrams corre- 
sponding to different topologies in the A;, — > pvr decays. In doing so, we give the expressions only 
for a certain representative set of diagrams in each class, with the rest following from appropriate 
substitutions. 



1. Factorization formulae for the color allowed emission diagrams 



For the first diagram in Fig. [5] (labeled as figure T\ ), which is a factorizable diagram and 
included only in the conventional pQCD approach, we have: 

J2 



Ti 



Gf 



^g^/Ajp / [dx] I [dx'] I dy [a,(t^i)]^VA,(x) 



,1 



,1 



1 



16Mlj(-Ci + C2)VubV:a + {-Cs + C4 + -Cg + Cio)VM]{-2x2 + (1 - 2x2)x', - x'^)^Uy) 



-32moMl{-C5 + Ce + -C7 + Cs)VM{4xi + 4x3 - x'3 - 3)</.f,(y) 



i/jpix'^ 



+ 



16Mlj(-Ci + C2)VM + (3C3 + C4 + 3C9 + c,o)vM]a + 4)^Uy) 
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+ ^16Mlj(-Ci + C2)K6Kd + (3C3 + C4 + -C9 + Cio)yffe-t^;j(2(xi + xs))<pUy) 
^ J b[db[J b2db2 J bsdbs J dOi j de2 exp[-S'^'ix,x',b, b')]Ko{Vm\b3\) 



167r2 
1 dzidz2 



zi{l-zi) \ \Z. 



X. 



(Bl) 



where the auxihary functions in the above expression are defined as 



A'^ = (1 - x'i)M? , = (X2 + x'^- X2x'^)Ml,C^' = X2x'2Ml, D^^ = xsx'^Ml, 



Z^' = A^^l - Z2) + 



Z2 



Zl(l - Zi) 



2^2 



(B2) 



Similarly, the factorization formula for the form factor /2 contributed by T\ can be written as 



Gp 



vr 



/aJp / {dx\ / [dx'] / dy{a,{f^)f^K,{x) 



,1 



1 



+ 



+ 



16Ml^ [(^C*! + C2)KfeKd + (^C-s + C4 + ^Cg + Cio)T46F;J(-2x2 + (1 - 2x2)x; - x',:)c\>i,{y) 



,1 



1 



-32moMl^(-C5 + Ce + -C7 + C78)yffeF;rf(4xi + 4x3 - x'3 - 3)</.i^(y) 



^,^(x') 



16Mlj(ici + C2)KfeKd + (^Cs + C4 + ^Cg + Cx^W:^\{2{x^ + X3))0^,(y) 



,1 



1 



-32moMl^(iC75 + Cg + ^Cy + C8)yt6^;rf(2(xi + X3))<Al/(y) 



Vj(x') 



^ 1 ^'idft'il 62<i62 1 &3*3 1 dd^ I de2exp[-S'^'{x,x',b,b')]Ko{Vm\b3\) 



167r2 

/■^ dzidz2 



X. 



Ti 



ZlO^^ ^|^l(\/^^)0(^2^') + + ^iVl(\/^^)]0(-^2^^)}- 

(B3) 

For the 25th diagram in Fig. [J] (labeled as T25), which is a non-factorizable diagram, we have: 



30 




FIG. 4: External W emission (T) diagrams for the A;, pn decay to the lowest order in the pQCD approach 
where the dots denote the weak interactions vertices. The two hard gluons are needed to transfer the large 
momentum to the light quarks in the initial state so that these two light quarks are coUinear in the final 
state. These diagrams are called Ti,T2, ...,736 in the text. 



16Mlj(-Ci - 2C2)VubV:a + {-Cs - 2C74 + -Cg - 2Cio)VM]{x2 - y){-x; + x'^ - y + l)0^,(y) 



(y)) 



+ 



+ 



+16moMli-C5 - 2C6 + -Cj - 2Cs)VtbV;,x'2{l - xs - y)Kiiy) 

16MlJ(^Ci - 2C2)VM + ilc3 - 2C4 + ^Cg - 2Cio)VtbV;a]{x2 - y){l - X3 - x'3 - y)<Pi{y) 



+16moMl(-C5 - 2Ce + -C7 - 2C8)VMx'2{l - xs - (y) 



(y)) 



i^^ix') 



32MlJ(-Ci - 2C2)VubV:^ + (-C73 - 2C4 + -Cg - 2Ci^)VtbV:^]{x2 - y){l - X3 - y)K{y) 



+32moMf (-C5 - 2C6 + -C7 - 2C8)VMx'^{y - X2)4>li{y) 



J b2db2 J bsdbs J bqdbq j dOi J d92exp[-S'^^^x,x' ,y,b,b' ,bq)] 



dzidz2dz3 



Zi)z2(l - Z2) 



X. 



T25 



-I 
2' 



725 I '7T25 I 



^725 I yT25 I 



7T25 



(B4) 



where the auxiliary functions in the expressions above are defined as 



^725 _ x'^(^X^ 

X^'" = {b2 - (63 - bq)z2 + bqZiil - Z2)f + | 



b„ + b„ziY + 



At 

2 , 21(1 - 2^1)22(1 - Z2) 2 



Z2ZZ 
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Zf^ = C{1 - Z3) + 



2^3 



^2(1 - Z2 



-[D{l-Z2) + 



Z2 



Zl{l - Zi) 



[A{l-zi) + Bzi]l 



=max{J\AT-^^^lJ\BT2^\,J\CT^^\,J\DT2'^\,oj,uj',ojg). 



(B5) 



Similarly, the factorization formula for the form factor /2 contributed by T25 can be written as 
f?' = GF:^^fAjpJ[dx] J[dx'] J dy[as{t^^^^)fi;Ad^) 

16Mlj(^Ci - 2C2)VM + ilcs - 2C4 + ^Cg - 2C,o)VM]{x2 - y){l - x^ - x'3 - y)^iAy) 



-16moM| r-C5 - 2Ce + -C7 - 2Cs)VtkVty2i'^ - X3 - y){cPfAy) " ^liv)) 



+ 



+ 



IQMl [(-Ci - 2C2)VubV:a + {-Cs - 2C4 + -Cg - 2Cio)VM]{x2 - y){-x; + x', - y + l)K{y) 



-l6moMU-C5 - 2Ce + -C7 - 2Cs)Vti>V:^x'^{l - X3 - y){(PUy) " 4>lM) 



32Mlj(-Ci - 2C2)KbKd + (3C3 - 2C4 + -Cg - 2Cio)VM]{x2 - y){l - X3 - y)Ki{y) 



-32moMU-C5 - 2Cq + -C7 - 2Cs)VMx'3{y - X2)(Pliiy) 



J b2db2 J bsdbs J bqdbq j dOi j d92exp[-S'^^'^{x,x' ,y,b,b' ,bq)] 



dzidz2dz3 



zi)z2{l - Z2) 



X. 



T25 



K,{^xJ-^zJ-^)e{zJ-^) + ^[M^xJ^^\zJ^^) + iN^i^xJ^^^\zJ^^\)]e{-zl-^ 



As can be seen from Eq. (jBSp . the color structures for the baryonic decays are quite different from 
that in the mesonic decays. Only the operators with the color indices the same as Oi can contribute 
to the non-factorizable emission diagrams in the two-body hadronic B meson decays. However, 
all the operators Oi{i = 1 — 10) contribute to the non-factorizable emission diagrams in the non- 
leptonic two-body bottom baryon Af, decays. In particular, the {V — A) ® {V + A) type operators 
have no effect on the non-factorizable emission diagrams for the hadronic B PP decays, if the 
emitted meson is a vr,?] or tj' . In contrast, both the (V — A) 1^ (V — A) and {V — A) f^i (V + A) 
operators contribute to the non-factorizable emission diagrams in their baryonic counterparts. 



(B6) 



2. Factorization formulae for the color suppressed emission diagrams 



For the first diagram in Fig. [5] (labeled as Ci), a factorizable diagram, we have: 



fi' = Gp^fUvjldx] J[dx'] j dy[as{t''^)f^A,{x) 

- 16Mf (-(C5 + Cr) + (Ce + CsWM K(x2(y - 2) + (y - l)y)«,(y) + (PIM) 



32 



-MA,(x2(y-l) + (y-2)y)<Af,(2/)] 



«(:r')-<(x')) 



32Ml (-(C5 + Cr) + (Ce + C^WM iv " 1) ["^0(^2 + y- l)(</'M(y) + </'M(y)) 



y" 62(i62y" fe^d^ay" bqdbqj dOiJ de2exp[-S^'{x,x',b,b')]Ko{VD^\b2 + b2-bq 



167r2 
^ dzidz2 



zi{l-zi) \ \Z[ 



7Cu 
'2 



2 - 



(B7) 



where the auxihary functions in the expressions above are defined as 

A'^ = (x'2 + y- x'2y)Ml.B^- = {x2 + y)Ml . C^^ = X2x'2M1 . D^^ = x^yMl. 



X. 



Ti 



A^^{l-Z2) + 
■- {b[ + Zib2f + 



Z2 



Zi{l - Zi) 
Zljl - Zi) 
Z2 



hi 



t^^ =max(v'|^^H,^/|5^H,^/|C'^M,^/l^''H,'^,^')• 



(B8) 



Similarly, the factorization formula for the form factor /2 contributed by C\ can be written as 

J2 



# = GF^fAjpJ[dx] J[dx'] j dy[as{f^)?i^K,{x) 

- 16Ml (-(C5 + Cj) + (Ce + C^))VM [mo{x2{y - 2) + (y - l)y){Ki{y) + <fii{y)) 



-MA,(x2(y-l) + (y-2)yX,(y)] 



{^'^{x')-^^{x')) 



32Mf (-(C5 + C^) + (Ce + Cs))^*^^;^ (2/ - 1) K(x2 + 2/ - miiiy) + <pM{y)) 



-MK,ct>Uy)] 



y 62^62 y 62^^62 y M^^y ^^^i y 'i^2exp[-5^H3;,x',6,6')]^o(^^^i&2+6'2-6g 



167r2 
1 dzidz2 





Zi(l - Zi)^ 



|^|Ki(yxF^)G(Z2^^) + l[J,{^|xf\zf\) + zN,{^lxf\zf\)]Q{-Z^^)Y 

(B9) 



For the 20th diagram in Fig. [5] (labeled C20 ), a non-factorizable diagram, we have: 



C20 



Gp 



27V3 



X <{ - 32M^^((C5 + C7) - {Ce + CsWM [MA.x^Kiiy) " 2mo(x3 - iX^(y) 
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-U 

-U 



3r 



— »- • < — 



5^ 



FIG. 5: Internal W emission (C) diagrams for the Ah pir decay to lowest order in the PQCD approach 
where the dots denote the weak interactions vertices. As in the preceding figure, the two hard gluons are 
essential to transfer the large momentum to the light quarks in the initial state. These diagrams are called 



j hidhi J bsdbs J bgdbgj dOi J d92eM-S'^''ix,x',y,b,b',bg)]Ko{VDC^\b2 + b'2-bg 



where the auxiliary functions in the above expression are defined as 



(BIO) 



Zf2o = A^^^'il- Z2) + 



Z2 



Zi{l - Zi 



■[B^'^''{l-zi)+C^^<'zi], 



Z2 



^1) 



fCao = max{\[\^^\, \[\B(^\, \[\C^l \[\D^\ 

, id, Oj', OOq). 



(Bll) 



Similarly, the factorization formula for the form factor /2 contributed by C20 can be written as 



C'20 



J2 ■ 



(B12) 
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3. Factorization formulae for the exchange diagrams 

For the 18th diagram in Fig. [6] (labeled as Eis), we have: 

^2 



El8 



Gf 



vr 



54^/3 



/aJp / [dx] / [dx'] / dyMt'^'^t^l^K^ix) 



i6moMl [(Ci - C2)VubV:^ + ((C3 + C9) - (C4 + Cio))x^tbF;d](y " mliiy) + 4>liiy)) 



+ 



16moMlj(Ci - C2)KfeK*d + ((C3 + Cg) - (Q + Cio))VM]{y - l)«,(y) + ^liiv)) 



2db2 J hdh J bgdbgj dOi j d02exp[-5^i«(x,x',y,6,6',6,)] 



167r2 



^ dzidz2 

Zl{l - Zl) 



X. 



\ ^1^1(^2''" ^2''")0(^2''") + 2 ['^l(V^2''"l^2''"l) +i^^l(\/^2''"l^2''"l)]0(-^2''") 

(B13) 



where the auxiliary functions in the expression above are defined as 



= (x'3 - l)Ml^, = (y - 1)(1 - x',)Ml,C^^^ = x'^iy - 1)mI, = x.yUl, 

ZElS ^ ^i?i8(l - Z2) + [B^^^l - Zl) + 



Z2 



max( 



(B14) 



Similarly, the factorization formula for the form factor /2 contributed from Ei^ can be written as 



jEia _ jEis 



(B15) 



For the 26th diagram in Fig. [6] (labeled as E2g), we have: 



ft 



Gf 



54^/3 



fAjp I [dx] I [dx'] I dy[as{t'''n]''i^A,{x) 



+ 
+ 



16moMlj(Ci - C2)VM + iiGs + G,) - (C4 + GwWM]il - y){<t>lM + <fiiM) 
16moMlj(Ci - G2)VubV:a + ((C3 + Cg) - (C4 + Cio))l^t6V^;d](l - y){<t^ii{y) + </>I/(y)) 
32Ml r(C75 + Gj) - (Ge + C8))ytbV;*,(y - 1)[Ma^{x[ - lX,(y) 



i;^{x') 
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FIG. 6: W exchange (E) diagrams for the A;, pn decay to lowest order in the pQCD approach where the 
dots denote the weak interactions vertices. As in the preceding figure, the two hard gluons are needed to 
transfer the large momentum to the light quarks in the initial state. These diarams are called Ei, E2, E^q. 



-moK-2)«,(y)+0i:,(y))] 



it^ I ^'i^^^'i/ ^'^"^^'^I ^^'^^''j "^^^ I d02eM-S''''{x,x',y,b,b',b,)] 



I Jo Zi{i - Zij 



1 dzidz2 



X. 



(BI6) 



where the auxiliary functions above are defined as 



A^^^ = {y- 1)(1 - x[)Ml,B^^^ = ix[ - l)Ml^,C^2« = 4{y - 1)m1^, D^^^ = xsyM^, 



[(6'2 + M+^l(^'l-^2-M] + 



2 , Zi) ^, ^ ^2 



Z2 



-{b[-b',-b,r, 
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= max ( J I I , J I 5^26 I J I I J I Z)S26 I 



(B17) 



f £^26 
/2 



Similarly, the factorization formula for the form factor /2 contributed by £'26 can be written as 
,2 



Gf 



TT 



54\/3 



fAjp / [fix] / [dx'] / [a,(t^2«)]VA,(x) 



+ 
+ 



i6moAflj(Ci - C2)VM + m + c,) - (C4 + c^oWt^vm " y)^^!^) + 



-moK-2)«,(y)+0^,(y))] 



{Ko{Vc^\b',\)9{C^'^) + '-l[Jo{Vc^\b',\) +iNoiVc^\b',m-C'''-)} /' ^£l^ 

2 Jo zi[i - zi) 



X 



E26 



4. Factorization formulae for the Bow-tie diagrams 



(B18) 



For the 17th diagram in Fig. [7] (labeled as -B17), we have: 



/i 



Gf 



fAjp / [dx] / [dx'] / [a,(t^^^)]2V'A,(x) 



- 16Mlj(-(C5 + C7) + {G, + C8))Fiby;rf](x'3 + x'i(y - I) - y + l)Kiy) 



i:Y{x') 



+ 



+ 



32moMf [(-Ci + C2)KfeKd + (-(C3 + C9) + (C4 + Gio))VM]{y - l){<PiM + fAM) 



-16Mlj(-(C75 + G,) + (Ce + G,))Vt,V:,]{x'^ + - 1) - y + l)</>^,(y) 



32M? [(-Ci + C2)K6Kd + (-(C3 + Cg) + (C4 + Gio))VM]x':i^ti{y) 



dzidz2dz3 



+32moM|j(-(C5 + Cr) + (Ge + C8))^t6V^ta(y - + </.I,(y))J 

^ / / .3^53 / / de, J dO. exp[-5-^^ (x, 6, b,)] [ - _ ^^^^^^^ _ 

(B19) 



Bl7| 

3 I 
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where the auxihary functions in the above expression are defined as 



1)M? , C^- = {y- 1)(1 - x[)Ml, I)^- = xax'gM, 



Bl7 



X. 



yBl7 
^3 



i{-b[ + h) - (-63 + bg)z2 - i-b[ + 63 - bg)ziil - Z2)f 



+ 



2:2(1 - Z2) 



A{\ - Z3) + 



((-63 + M-^i(-&'i + &3-M)' + 



zi(l - 21)2:2(1 - Z2) 



Z2. 



2:2(1 - 2:2) 



5(1-22) + 



Z2 



2:1(1 - Zi) 



Z2Z-i 

[C{l-zi) + Dzi] 



i-b[ + 63-6, 



(B20) 



Similarly, the factorization formula for the form factor /2 contributed by fig. Bn can be written 



as: 



J2 



Gf 



216^/3 



/aJp / [dx] / [dx'] / dy [as{t^'n?^^ASx) 



32moMf[(-Ci + C2)VM + ("(^3 + Cg) + (C4 + Cio))Vtby;,](y - l){(t>lM + fuiy)) 



+16Mlj(-(C75 + C7) + (Ce + C78))T4,y;,](4 + x',{y - I) - y + l)</.^(y) 
16Mi[(-(C5 + C7) + (Ce + C^))VM]{x'^ + ^(y _ 1) _ y + l)0^^(y) 



32M1J(-Ci + C2)VubV:, + (-(C3 + C9) + (C4 + C^o))VM]x'^i^Uy) 



-32moM|j(-(C75 + C7) + (Cg + C8))14,y,a(y - l)(./.f^(2/) + fuiv)) 
32^ / ^'^"^^'i/ I ^"'^^"j "^^^ I ^^2exp[-5^i^(x,x',y,6,6',6,)] 



V'J(x') 



(izidZ2(i23 



Bl7 



\\Z: 



Bl7| 



zi(l - 2:1)22(1 - 2:2) 

K^{^xf^)eizi^^) + |[Ji(Vx3^-|zf-|) + ziVi(v/xf^^^)]e(-Z3^-)}. 



(B21) 



For the 19th diagram in Fig. [7] (labeled as -B19 ), we have: 



/i 



Bi9 



+ 



^fAjpJ [dx] J [dx'] J [a,(t^i«)]VA.(x) 

- i6Mlj(-(C5 + C7) + (Q + C8))i4fey;rf]4'/'M(y) 

- 64moM^J(-Ci + C72)K6K"d + (-(C3 + Cg) + (C4 + Cio))VtbVt2<pUy) 



i^l{x') 



-16Mlj(-(C5 + Cj) + (Ce + C78))yifey,a40^(y)J V'^(x') 
32Mlj(-Ci + C2)VubV:a + (-(C3 + Cg) + (C4 + Cio))VM]x'2^ti{y) 
-64moM|j(-(C5 + C7) + (Ce + C8))Vtby,a(x'i - 2)</>^(y) 
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d- 










i 








FIG. 7: Bow-tie (B) diagrams for the pn decay to lowest order in the pQCD approach where the 

dots denote the weak interactions vertices. As in the preceding figures, the two hard gluons are needed to 
transfer the large momentum to the light quarks in the initial state. 



^ J b[db[J b'^db'^J bgdbgj d9i J d02exp[-5^«(rc,x',y,6,6',6,)] 

{Ko{Vc^\b,\)6{C''^^) + '-l[No{Vc^\bg\) + iKo{Vc^\b,m-C''^^)} f' -^^'^''^ 

^ Jo Zi[i 



X. 



Big 



(B22) 



where the auxiliary functions in the above expression are defined as: 



A' 



c[ - 1)M? ,C^« = x'.iy - 1)M? ,1)^" = xax'gM 



Z2 



Zl{l - Zi) 



-zi) +Z)^i«zi], 



Z2 



^1 



Bi9 
2 

= max ( a/ I I , J 1 5^19 1 , J I C^w I , J I L>^i9 1 , cj , , cj, 



(B23) 



Similarly, the factorization formula for the form factor /2 contributed by -Big can be written as: 

- UmoMi\{-Ci + C2)VubV:^ + (-(^3 + Cg) + (C4 + Ci^m^V^^^t^Uy) 



/f« = Gp^fKj:, l[dx] l[dx'] I dyiasit^'^nfMx 
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+16M? [(-(^5 + Cj) + (Ce + Cs))VM]x',Kiy) 



16M? [(-(Cs + Cj) + (Ce + CsWMW.^Uy) 



32M1J(-Ci + C2)K6Kd + (-(Cs + Cg) + (C4 + Cio))I4fey;,]x^V'f,(y) 



167r2 



^W^'i/ b2db'2j bgdbgj dOi j d92 exp[-S^^^ {x,x' ,y,b, b',bq)] 



{KoiVc^\bg\)e{c''^^) + '-^[No{Vc^\b,\) + tKoiVc^\bgm-c''^^)} 



dzidz2 
zi{l - Zl) 



X. 



Big 



Bl9\ 
2 I 



(B24) 



5. Factorization formulae for the penguin annihilation diagrams 



For the 14th diagram in Fig. [8] (labeled as P14), we have: 



Pl4 



54^/3' 



/aJp / [rfx] / [dx'] I dy [a,(t^")]2V;A,(x 



48MlJ(C3 + C4) - -(Cg + C,o)]VM{MA,iy - l)<Pij{y) + moicf>Uy) - 4>li{y))) 



1 



1 



1 



32Mf [(C5 - -C7) + -(Ce - -CsWMyiMAAKiy) + m,{<Pi,{y) - 0^/(2/))) 



2 

-32Ml \\{C^ - \cj) + (Ce - ic8)]ytfey;rf(MA,x^0^(y) + mo(y - 2)</.r^(y) + mo#If(y)) 



<(x') 



+ 



1 



48M1J(C3 + C4) - -(Cg + C^,)]Vt,V:^{M^,{y - l)0^,(y) - mo(0i:,(y) - 
+32Mlj(C5 - ^Cr) + \{C^ - \c^)]VMy{M^,x'2(t>Uy) " "^o(<Alf (y) - (y))) 



.1 



1 



1 



+ 



+32M^J-(C5 - -C7) + (Cg - 2C8)]^t6V^;d(MA,x^0^,(y) - m^{y - 2)<p'iiM - moy<PiM) 

96moM|j(C3 + c^) - \{c^ + Ci^)]vMy{4>M{y) - <Pliiy)) 
+64Mly[{c, - \c7) + \{c^ - \c,)]vM{y " ^)'PUy) 



^ y b2db2 j bgdbgj b'2db'2j dOi j de2exp[-S^'^{x,x',y,b,b',bg)] 
{Ko(Vc^|62|)e(C^") + f [iVo(V^C^|62|) + ^i^o(V^C^|fo2|)]^(-C^-)} ^^^^ _ 



dzidz2 
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(B25) 



where the auxihary functions in the above expression are defined as: 

= x'^MI^,bP-^ = (1 - x'^y)Ml^,CP^^ = X2x',Ml, D^^^ = x[iy - 1)MI 



Z2 



XP' = [(62 + b'2) - z,bf + ML^5, 2 



Zl{l - Zl) 
2 



^2 



q ' 



(B26) 



Similarly, the factorization formula for the form factor /2 contributed by P14 can be written as: 
/f^^ = GF^hJpJidx] j[dx'] I dy [asit^'nf i^A,ix) 



1 „ , 1 1 



+32M1J(C5 - -Cr) + 2 (Ce - -C8)]Fi6F;,y(MA,x^0^,(y) + mo(<Ai^(y) " Kiiy))) 
-32Mlji(C5 - ^Ct) + (Cg - ^C8)]ytfcy;rf(MA,x^</.^(y) + mo(y - 2)<A^(y) + mo#If(y)) 
48Mlj(C3 + C4) - ^(Cg + C^o)]VM{Mj,^{y - l)<pUy) - mo«f(y) " ^liiv))) 



-32Mlj(C5 - ^Cr) + ^(Ce - ^C8)]yiby;,y(MA,x'2</.^(y) - mo(<Al^(y) " 0M(y))) 



+ 



+64M5 y[(C5 - ICr) + ^(Cg - ^C8)]I4bF;,(y - l)<A^(y) 



ib^Lv-^ 2 ~ " ■ 2'~" 2 
-1 ,„ 1 „ . 1 



-64M1J-(C75 - -Cr) + (Cg - -Cs)]VM{x', - l)</)^(y) 



{Ko(^^|62|)e(C^") + Y^No{^/\C^\\b2\) + ii^o(V^C^|62|)]e(-C^")} 



dziC?Z2 
Zl{l - Zl) 



X. 



Pl4 



l^j A'i(v/Af"zf»)e(zf") + ^iJil^xf-'izf"!) + ijv,(\/xf"|zf"|)]e(-z|'" 

2 



(B27) 



For the 26th diagram in Fig. [8] (labeled as P26), we have: 



P26 



Gp 



IT 



27V3' 



IaJp / [dx] / [dx'] / [a,(t^^«)]VA,(x) 
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b. 

U- 




FIG. 8: Penguin annihilation (P) diagrams for the Af, — > pir decay to lowest order in the pQCD approach 
where the dots denote the weak interactions vertices. As before, the two hard gluons are essential to transfer 
the large momentum to the light quarks in the initial state. These diagrams are called Pi, P2, ^36- 



16Mlj(C5-^C7 



1 



1 



1 



+ 



2 

1 ^ . 1 



-16MlJ(C5 - -C7) - (Ce - -Cs)]VtbV;a{MA^x[<Ptiiy) + 2mocPM) 



16M4J(C5 - -C^) - (Ce - -Cs)]VtbV:^{Mp,^x\<l,Uy) + 2mo{x', - 2)0^,(y)) 



i;^{x') 



^ J b[db[J b2dh2 J b,db,J dOi J d02exp[-S''^'{x,x',y,b,b',bg)] 

2 " V Jo - zi) 



dzidz2 



(B28) 



where the auxiliary functions in the above expression are defined as: 

= -x'lM? = {x'^ - l)Ml, C^^s = X2x'^Ml,D''^^ = x[{y - 1)MI 



Z2 



Zi(l - Zi) 



[B^^^'il- zi) + C^^^zi] 



xf- = [(&; + b,) - z,b2f + 



Z2 
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i^26 = max{J\AP^^>\,J\BP^6lJ\CP^^\,J\DP^^\,uj,u)',ojq). 



(B29) 



fP26 
J2 



Similarly, the factorization formula for the form factor /2 contributed by P26 can be written as: 
^2 



vr 



GF^fAjpjidx] j[dx'] jdy [a,(t^2«)]VA,(x) 

- 16Mlj(C5 - ^Ct) - (Ce - \Cs)]Vt,V:^{M^^x'^<pi,{y) - 2mo{x', - 2X^(y)) 



1 



1 



+ 



16Mlj(C5 - ^Cr) - {C, - ^Cs)]VMiM^,x[^Uy) + 2mo(4 - 2),/.^,(y)) 
+16Mlj(C5 - ^^7) - (Cg - ^Cs)]VthV,mA,x[c^Uy) - 2mo<,(y)) 



^ 1 ft'ldfo;! 62(i62 1 bgdbgj dOi J d02exp[-5^2«(x,x',y,6,6',6,)] 



dzidz2 



P26 



(B30) 



6. Factorization formulae for the three-gluon-vertex diagrams 



Now, we can focus on the hard amplitudes contributed by the topological diagrams shown in Fig. 
[3] with the insertion of the three-gluon-vertex, which have been grouped in Fig. [9l It needs to be 
pointed out that the insertion of the three-gluon-vertex to the external and internal W emission di- 
agrams, namely the diagrams GTi and GCi {i = 1 — 4) in Fig. [9]have null effect on the decay ampli- 
tude, since the color factors in these diagrams are proportional to eijkf-i' j'k' f"'^'^{T'^)i' i{T^) j' j{T'^)k'k^ 
which equals zero taking into account the symmetry property of the structure constant This is 
also the reason why the Feynman diagrams with the three-gluon-vertex are neg lected in computing 
the hard amplitudes for the semi-leptonic decays of the A;, baryon 
For the 1st diagram in Fig. [9] (labeled as GEl), we have: 
^2 



15 



13, 



39| 



GEl 



Gf 



vr 



24^/3' 



/aJp / [dx] / [dx'] / dy [a,(t^^i)]VA,(x) 



16Mf ((Cs + Gj) - (Ge + GsWMiM^.iSx'^ + X3)cPi,{y) 



+mo((3 - 2x3)<^(2/) - (1 - 2y - 24))0lKy) 



V'J(x') 
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d u 




(GPl) (GP2) (GPS) (GPi) 



FIG. 9: Feynman diagrams responsible for the Ah — s- pi: decay with three-gluon-vertex to the lowest order 
in the pQCD approach where the dots denote the weak interactions vertices. 



2 " " Jo Zi(L 



dzidz] 



2 



zi{l - Zi) 



where the auxiliary functions in the expression above are defined as 

Zi{l - Zl) 

X^'^' = [{-h, + h', + h,) - z,h,f + 



i^^i = max{^\AGJ^^,^\BGEil^\cGEi\^^\jjGEi\^^^J^^^Y (B32) 

Similarly, the factorization formula for the form factor /2 contributed by GE\ are written as: 

jGEl ^ _fGEl^ 



2 



(B31) 
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For the 3rd diagram in Fig. [9] (labeled as GE3), we have: 



GE3 



-Gf 



24^/3 



fAjp I [dx] I [dx'] I dy [a,(t«^i)]VA.(^) 



8Ml[(Ci - C2)VubV:d + ((^3 + Cg) - (C4 + Cio))VM] 



+ 



x(2AfA,40M(y) + "io(30i;(y) + (1 - 2y)<Pi,{y)) 

-SMliiC + Cj) - (Ce + C^))VM{2M^,x'2<t>Uy) + ^oi^^niv) + (1 " ^yW^M) 

- 8Mlj(Ci - C2)F„feKd + ((^3 + Cg) - (C4 + Cio))l^tfey;,] 

x(2Ma,x^0^,(2/) + mo(3<,(y) + (1 - 2y)(/.I,(2/)) 

+8M|^((C5 + C7) - (Ce + Cs))VM{2Mp,^x'^4>i,{y) + mo(3</.:^(y) + (1 - 2y)(g,{y)) 
16Ml((C5 + Cr) - (Ce + Cs))^^,,^;,!^^ - l)(MA,(2y - l)<t^Uy) 



-mo(3(y - iWhiiy) + (1 + y)0M(y)) 



^ 1 6'2d6'2y' 63^63 1 hqdbqj d6i J d92exp[-S^''^{x,x',y,b,b',bg)] 



^ dzidz2 

Zi(l-Zi) 



(B34) 



where the auxiliary functions in the above expression are defined as: 



A 



GE3 



ix',-l)Ml,B 



GE3 



GE3 



x'2{y-l)Ml,D 



GE3 



X3yMl 



ZGE3 ^ j^GE3^^ _ + 



Z2 



Zi{l - Zl) 



GE3. 



X<f^' = [{b2' + b,)-Z,b,]'+''^^-''^-' 



'■2 
J.GE3 



Z2 



t^^' = max{J\AOE3\J\BGE3\J\cGE3\J\DOE3\^^^^'^^^~) 



(B35) 



Similarly, the factorization formula for the form factor /2 contributed by GE3 can be written 



as 



GE3 



-Gf 



TT 



24^/3 



/aJp / [dx] I [dx'] I dy[as{t''^^)]^^PA,{x) 



8MU{Ci - C2)VubV:, + m + Cg) - (C4 + CwWtbVth 



x{2MKAnAy) + moi^KM + (1 - 2yWM{y)) 

-SMUiC^ + Gj) - (Ce + G^))VM{2MA^x'^4>i,{y) + mo(3</.:P (y) + (1 - 2y)</.^(y)) 
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+ 



8Mlj(Ci - C2)VubV:, + ((C3 + Cg) - (C74 + C^oWM] 



x(2MA,40M(y) + mo{'i(l}%i{y) + (1 - 2y)(/<'l,(y)) 



T 



+8M1J(C5 + C7) - (Ce + Cs))VMm^^x'^<l)Uy) + "^oCSc^-IfCy) + (1 - 2y)(/<i,(y)) 
- 16Ml f(C5 + Cr) - (Ce + C78))14bF;,(x;i - l)(MA,(2y - 



-mo(3(y - l)(t>li{y) + (1 + y)0M(y)) 



^ j b'^db'^J badbs J bgdbgj dOi J de2eM-S^''Hx,x',y,b,b',bg)] 

{KoiVJGE^lb'mA''''') + ^[Jo(V^4^|6'2l) +iiVo(vt4^|6'2l)]^(-^''''')} , 

2 " " Jo zi(L - zi) 



dzidz2 



GE3\ 
2 ) 
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